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Abstract 

I develop the simplest geometric-discretized analogue of two dimensional scalar field 
theory, which qualitatively reproduces the trace anomaly of the continuous theory. The 
discrete analogue provides an interpretation of the trace anomaly in terms of a non- 
trivial transformation of electric-magnetic duality-invariant modes of resistor networks 
that accommodate both electric and magnetic charge currents. 
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1 Introduction 



Finding new formulations of familiar concepts frequently proves beneficial in physics. This 
may be motivated by computational needs, as in lattice field theory, or by a hope to gain 
new fundamental insights. Since the 1960s the idea of geometric discretization has lured 
physicists with promises of new ways of looking at field theory and quantum gravity. Af- 
ter Regge calculus PQ, a discretized theory of gravity, a discrete understanding of Abelian 
Chern-Simons theory [2], BF field theory [3] and chiral fermions [4| has been attained. Not 
surprisingly, geometric discretization has been most successful in modelling topological con- 
cepts and theories. Meanwhile, progress in discretizing non-topological objects has been less 
convincing. Although Desbrun et al. [5] provided explicit discrete analogues of geometric 
constructs including the Hodge star and the Lie derivative, the complexity of their formal- 
ism makes it difficult to derive insights from it. For example, at present there is no discrete 
derivation of the well-known trace anomaly in two dimensions. 

Here I develop a new, simple discretization of scalar field theory and use it to derive 
the form of the trace anomaly from first principles. I treat the discrete theory as defined 
independently of its continuous counterpart and work only with constructs which are intrinsic 
to the discrete theory on a triangulation. The starting point of the discretization is the 
identification of p-forms on a manifold with p-cochains on a triangulation, and thereafter the 
calculation proceeds along a path which is essentially selected by consistency requirements. 
As such, the present work exhibits the simplest discrete analogue of two-dimensional scalar 
field theory on a metric space. 

At the end of the calculation a new way of thinking about the trace anomaly emerges. 
The manifold (triangulation) is viewed as an electric-magnetic network and conformal trans- 
formations are found to act non-trivially on the relevant charge densities. Appealing to the 
correspondence with the continuous case, the final form of the discrete trace anomaly may be 
viewed as singling out a preferred, smeared out notion of curvature on a finite triangulation 
of a two-dimensional manifold!]] More on this in the discussion, which follows the review of 
Sec. [2] and the calculations of Sec. El 



2 Review of the continuous theory 

I begin with a brief review of the scalar theory on a smooth compact orientable two- 
dimensional manifold M without boundary, based on [7]. This will help to fix all relevant 
conventions and provide a necessary reference point for setting up a dictionary between the 
continuous and discrete constructs. 

A manifold M is endowed with the de Rham complex Q*(M) 

^ n°(M) ^ Q}(M) % n 2 (M) % 0, (1) 

where the maps d r satisfy d r+ id r = 0, and the wedge product 

A r : VL r (M) x n 2 -' r (M) -> C 

7 (2 - r) ) = | M aWA r 7 (2 - r) = H) r / M 7 (2 ~ r) A 2 - r a (r) . { ) 

1 For an alternative treatment, see Ko and Rocek [BJ. 
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The introduction of a metric g^ v induces the Hodge star, a natural isomorphism between 
Q r (M) and Q 2 ~ r (M), which in turn induces a canonical inner product on Q r (M): 

* r : VL r (M) -> VL 2 ' r (M) 

where in the second line I suppressed the superscripts indicating the degrees of the forms. I 
take the manifold M to be Riemannian, so the Hodge star satisfies: 

*2-r 0* r = {-iy. (4) 

The Hodge star allows one to introduce operators 8 r+ % adjoint to the exterior derivatives d r : 

5 r+1 : Q r+1 (M) -> fT(M) 

8r+l = *2-r ° dl~r ° *r+l • 

The resulting structures, the de Rham complex together with its adjoint, are summarized in 
the diagram below: 

i n°(M) $t fi^M) i ^ 2 (M) i 0. (6) 

#o 5i 52 i 

A 0-form field $ coupled to a gravitational background on a 2-dimensional manifold M 
without boundary has the action: 



S = l [ d & A 1 M $ = ~<d <I>, . 



(7) 



The appearance of *i in (GO) breeds a dependence on g^ u and is the ultimate reason for the 
theory not being topological. The effective action 



W[g] = -log j[d$]e 

suffers from an anomaly: 



-S[*a] 



bg^ u {x) 2 
The small time diagonal heat kernel takes the form 

\\m(x\e^\x) = -L + JLfl(aj) + 0(e) , (i ) 
e-^o 47re 247r 

which after the addition of a local counterterm 5 — > S — A/8ire yields 
The anomaly, integrated over M, gives: 

/ ^-jj = -7.X- (12) 
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3 The discrete theory 

I now set out to find the analogues of eqs. (171 fT2l) for a theory defined on a finite triangulation Q 
of the manifold M. Q is most easily thought of as a graph G(V, E) containing \V\ vertices and 
\E\ edges, embedded in M in such a way that no two edges cross. For technical convenience 
I shall assume that G does not contain loops and the number of edges E is even. The 
embedding partitions M into \F\ faces, which are (homeomorphic to) triangles when G 
corresponds to a triangulation but may be (homeomorphic to) arbitrary polygons for a more 
general G. The number of faces is such that the Euler characteristics agree: 

Xg = \V\ + \F\-\E\= X m. (13) 
It is natural to associate p-forms on M with p-cochains of Q. Because we are interested in 



modelling a real scalar $, I shall take Q° xeM — > Q®ev = ^ 


I. Thus: 




Q°(M) < — ► R v 






n\M) < — ► r e 




(14) 


tt 2 (M) < — >R F 







This set-up suggests a canonical correspondence between the exterior derivative on M and 
the co-boundary operator on Q. To define the latter, assign to each edge e G E an arbitrary 
orientation, that is to say, think of one of the two vertices incident to e as "initial" and of 
the other one as "final" . Then the \E\ x \V\ directed incidence matrix M G given by 

{—1 if v = init(e) 
1 if v = fin(e) (15) 
otherwise 

acts on ~R V as a coboundary operator [8] on the space of 0-cochains: 

M G : R v -+R E < ► d : Q°{M) -> tt^M) (16) 

The image of M G in R E is called the bond space of G, denoted 13(G). If we think of G as 
specifying an electrical network, 13(G) is the space of all possible potential differences. The 
orthogonal complement of 13(G) in ~R E is the cycle space C(G), the space of all circular flows 
in G. The analogy with an electrical network will be useful to us later. 

If G represents a triangulation, the correspondence f|T6|) is canonically extended to d\ as 
follows. Denote every edge e G E by [ij], where Vi = init(e) and Vj = fin(e). In a similar way, 
denote every face / G F by [ijk], where Vi, Vj, Vk are incident to / and their ordering agrees 
with the orientation of the underlying manifold M. Permutations act on [ij], [ijk] like parity, 
e.g. [ij] = —[ji] and [ijk] = [jki] = —[jik]. We define the linear operator N G : M. E — > R F 
via: 

(N G uj)([ijk]) = u([ij]) + uj([jk]) + u([ki]) for u G R E (17) 
Then by construction N G M G = 0, which mimics did® = under the correspondence: 

N G : ~R E — > R F < — ► d x : Q\M) -> Q 2 (M) . (18) 
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It is easy to see that the image of N G spans all of M F except the constant mode. This agrees 
with the Hodge decomposition theorem, which stipulates that 



dimfi 2 (M) = dimdifi^M) + b 2 . (19) 

The only remaining non-trivial consistency condition is the antisymmetry of the wedge prod- 
uct A , which is made possible by the evenness of \E\. Consider any matrix of wedge products 
A 1 (G): 

A 1 ^)^ = K u 2 ) (20) 
Then the analogous A°(G) is determined by enforcing Stokes' theorem: 

/ cW 0) A 1 7 (1) = - / « (0) A dn (l) — (M G ) T A 1 (G) = -A°(G) N G . (21) 

JM JM 

The above definition is supplemented by requiring that the matrix element of A°(G) evaluated 
in the constant modes M F \Im(N G ) and ker(M G ) be non-zero, so as to ensure non- vanishing 
area of the manifold. In summary, the topological data of the manifold M find well-defined 
analogues in the triangulation Q . 

However, the introduction of a metric interrupts this complacent state of affairs. The 
isomorphism ([3]) requires Q r (M) and Q 2 ~ r (M) to be of the same dimension. Under the 
correspondence (fT4"j) this translates into |V| = \F\ and the only triangulation which satisfies 



this is the tetrahedron triangulating a sphere. If we relax the condition of working with a 
triangulation, the map ffl8|) will be lost; meanwhile, \V\ = \F\ remains a cumbersome and 
non-trivial condition^ A more general possibility is to replace G with the union of the graph 
and its dual, GUG, with the understanding that the resulting p-cochains should correspond 
to a reduplicated space Q P (M) © Q P (M) [2J. In this paper I follow the latter strategy but 
abstain from the onerous detail. 

Consider two copies of the scalar field theory defined in (J7J). The fundamental scalars 
$! and J?2 are naturally associated with the 0-cochains of G U G whose space is isomorphic 
to M^®^. The previous notion of the exterior derivative d <-» M G naturally extends to G. 
First one identifies the sets E and E, that is, M. E 3 e = e G and M. E = M. E . Then one 
extends the orientation on G to G in the following fashion: for an edge e G G separating 
faces fi, fj G F we say that init(e) = fi if in crossing from to fj in G one passes fin(e) on 
the right and init(e) on the left; otherwise fi = fin(e). This prescription is well-defined by 
virtue of the orientability of the manifold M in which G is embedded. Then one forms an 
\E\ x \F\ (i.e. \E\ x \F\) directed incidence matrix M G as in (1151) : 




if / = init(e) 
if /=fin(. 
otherwise 



M° f ={ 1 if/ = fm(e) (22) 



The matrices M G and M G can be concatenated to form the collective incidence matrix 
M GuG , whose dimensions are |i£|x(|V| + |.F|). It is easy to see that M G sends M. F into 
C(G), that is, B(G) _L 13(G). In other words, a potential difference in the network G is 



2 This program was pursued in [S]. 
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necessarily a circular flow in the network G. However, the converse is not true for non-zero 
genera of M. There are 2g linearly independent flows which are circular in both G and G, 
i.e. dimC(G) D C(G) = 2g. This observation will be of essential importance later. 

Because we are working in G U G, the triangular character of the underlying graph is in 
general lost and we do not have the luxury of extending ( fTTI) to the present case. However, 
since the action for two scalar fields reads 

If If 1 1 

Si2 = - #iA*#i + - / d$ 2 A*d$ 2 = d$i) + -(d$ 2 , d$ 2 ), (23) 

2 J M 2 J M 1 1 

finding an explicit operator representing d\ is not necessary. Similarly, an explicit repre- 
sentation of A 1 is redundant as the only combination that shows up in the action is A 1 *!. 
Thus, in the following I shall not be concerned with searching for discrete analogues of all 
the standard machinery of algebraic topology and geometry. Instead, I shall limit myself to 
making a simple ansatz for the discrete representation of f M . . . A 1 *\ ... = (...,.. .) This, 

complemented with a notion of area and the correspondence M GuG d , will be sufficient 
for determining the trace anomaly. 

The introduction of a metric is in order. This is done by assigning to each edge e G E a 
length g e . I then take the inner product on M. E to be: 

J - (...,...}« 

Jee' = 9eS ee ' (24) 

To complete the mapping, a notion of area is necessary. This is defined as: 

^ M e&E v€V 



where the degree of each vertex d v is the sum of the lengths of the edges incident to v. The 
correspondence ([25]) has been considered in mathematical literature (viz. Sec. 3.4 in [10] ). 
although it is not unique. Because differences between the present approach and others (for 
example, [6]) may ultimately be traced to the choice of metric, eqs. ( 12"4H2"oT) represent a key 
step in the present development. I shall return to this point in the discussion. 

The correspondence (I24ti25|) is well-motivated. From an abstract point of view, no prior 
restrictions constrain the discrete analogue of the Hodge star so the initial choice of mapping 
A <-> vol(Cr) is unrestrained. After that, the inner product (. . . , • • is expected to have 
the same scaling as (. . . , . . .)(°) and exhibit an appropriate metric-related notion of locality 
in IR 5 . This essentially determines (...,.. .)^ in a unique way. 

On a more intuitive level, the analogy with an electrical network is helpful. The zero 
modes in M. v and M F , the respective constants, can be thought of as specifying constant 
densities of excess $i and $2 charge in the network. If we envision the charges as residing 
along the wires of the network, then it is a canonical choice to demand: 

Qtotal = CO A <-> A 2J 9e = Qtotal , (26) 

eG-B 

where <7o and Ao denote the respective constant charge densities. The vectors in M E specify 
flows in the network. Since (1, oc Qtotal, it is natural to require that (l e , le)^ be 
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proportional to the amount charge flowing in e G E, which in turn is proportional to the 
length of the edge g e . This motivates (1241) . 

With the definitions of M GuG do and (. . . , • • in place, the action fl23|) reads: 

S l2 = ^ T (M GuG ) T J M Gu % , (27) 

where \l/ G M^®^ represents both $x and $2 and lives in the direct sum of their respective 
field spaces. The effective action is given by 



W 12 [g] = -log J[d%e- Sw . (2t 



There are two zero modes which need to be regulated, one for each of $1, $2- They 
are given by the constant modes on V and F, respectively. Using the electrical networks 
analogy, it is natural to pick the following regularization scheme. Suppose the network is 
capable of sustaining static charge densities only up to a maximal value X l max , i = 0, 1. 
In terms of the normalized zero mode ip l = (vol(G)) -1 / 2 , the charge density is given by 
< Aq = CqV'o — c o( vo U^)) ~ 1 ^ 2 — Knax- Then the integral over each zero mode returns the 
value of (c ) max = X max {vo\{G)fl 2 . Thus, eq. (ESD yields 

W 12 [g] = -logvol(G) + logdet'{(M GuG ) T JM GuG ] (29) 

up to constant terms involving \ l max . 

The second term in ( 12"9"|) cannot be split into logarithms of determinants of the factor 
matrices because the rank of M GuG is lower than that of J. The obstruction to doing so 
lies, therefore, in the zero modes of (NL GuG ) T : M> E —> (My eF ) v . I proceed to regularize these 
modes. 

The strategy will be analogous to the regularization of the zero modes of 'M. GuG . The 
relevant zero modes are flows which are circular both with respect to G and G. There are 
\E\ — (\V\ + \F\ — 2) = 2 — x linearly independent such flows. Suppose that for a particular 
zero mode Ti G C(G) flC(G) the underlying electrical network can only support charge fluxes 
not exceeding a maximal value Z™ 11 '. This means that the flux density over an edge e is 
bounded by: 

< CiT f (e) (5>(e)V)~ < XT X n(e) (30) 

eg_E 

1 /2 

and the volume of integration over q is ( ^2 eeE Ti(e) 2 g e ) T™ ax . The regularization allows 
us to treat the determinant in (12TJ1) as if it were non-singular: 

det'{(M GuG ) T JM GuG ] = detJ det |m GuG (M GuG ) t } . (31) 

The second determinant is not primed as it includes the regularized zero mode contributions. 
We have: 

det |M GUG (M GUG ) T } = det' |m GuG (M GuG ) t } jf Z[ nax ( ^ r^efg^ ^ (32) 

i=l e&E 
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The index % runs over an orthonormal basis T{ of C(G) D C(G). 

Although the first factor is independent of the metric and will merely contribute a con- 
stant term to it is amusing to note that it is blessed with a beautiful combinatorial 
interpretation. If neither G nor G contains loops, then the determinant factorizes over (jt, G 
and according to the matrix tree theorem [TT] 

det ' {M G (M G ) T } = \E\ #{ spanning trees of G} (33) 

with the analogous equality holding for G. This result, first demonstrated by Kirchhoff in 
1847, is useful in calculating total resistances in electric networks which are neither parallel 
nor in series. 

Putting all the results together and dropping constant terms, the effective action takes 
the form 

W[g] = -\ log ( £ 9e) + \ E log ( E T ^?9e) + \ E lQ S 9e • (34) 

e&E i=l eG-E eG-E 

Here I have restored the factor of 1/2 coming from the fact that Wi% contains contributions 
from two scalars $1, $2- The conformal variation of W is given by 



9, 



5W[g] 1 g e 1 ^ Ti (e) 2 g e , 1 



e 



1 Tile) q e 1 ,„^ s 



As in the continuous case, the last term may be cancelled by introducing a local ^-independent 
counterterm S — > S— 1/2 Tr log J. The first two terms, however, form a topological invariant: 

2l9e^ = - A X. (36) 

The correspondence with the discrete result (flT| [T2|) is imperfect - the coefficients do not 
agree. Nevertheless, eq. (135]) singles out a preferred discretized notion of curvature, distinct 
from that of [T]: 



A 9e (lg^|L- 1 ). ,37) 



4 Discussion 



This paper presents the shortest route to discussing the trace anomaly in discretized settings. 
As a starting point, the canonical correspondence between p-forms and p-cochains on a 
triangulation was adopted, along with the natural definition of the th exterior derivative. 
After that, the metric was introduced via the Hodge star, leading to a doubling of the field 
space for consistency. Every construct introduced in this development is strictly necessary 
for a discussion of the trace anomaly, and conversely, the construction is free of redundancies. 
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It is in this sense that the presented theory is the simplest possible in the equivalence class 
of consistent discrete two-dimensional scalar field theories exhibiting trace anomalies. 

The derivation involved one arbitrary choice - that of the discrete Hodge star *x (and 
therefore the metric) in eqs. (I24H25I) . The correspondence adopted in this paper was mo- 
tivated by an effort to treat the objects intrinsic to the triangulation - vertices, edges and 
edge lengths - as fundamental, and avoid making recourse to an explicit embedding E It 
should be borne in mind, however, that other equally motivated choices are possible. Of 
these, a notable one is that employed by Ko and Rocek in a related paper [6] , which follows 
the tenets of Regge calculus. In it, the discrete theory is implicitly assumed to live on a 
piecewise linear space consisting of a union of flat triangles, so that the area of the manifold 
does not conform with eq. (125]) . This leads to an effective action, whose conformal varia- 
tion features the conical curvature and which differs from eq. ( 1341) . However, the general 
approach presented in this work is robust. With a different input metric to replace eqs. (1241 - 
[2"5l) . it should in principle reproduce the result of [B|. The numeric disagreement between the 
discrete anomaly (IBljj) and its continuous counterpart (1T2"|) is likewise easily fixed by a trivial 
redefinition of the metric g e — > g e 3 ^ 2 in eq. ( 1241 . It would be interesting to understand the 
exponent —3/2 in detail, as it provides a quantitative bridge between two-dimensional scalar 
theory and electric-magnetic resistor networks. 

The latter analogy is a nice bonus of the calculation presented in Sec. |3j We interpret 
the zero modes of M, as charge densities and the zero modes of M T , Tj's as circular flows 
in a network of resistors. The network is electro-magnetic as both electric and magnetic 
charges propagate in it. This is reflected in the fact that the field \1/ has two independent 
zero modes, one living on G and the other on G. Indeed, upon the exchange G ^ G we 
have: _ 

B(G) C(G) m) 
B(G) ^ C(G). y ' 

This operation is easily seen as an interchange of the two-dimensional divergence and curl, 
which effects an electric-magnetic duality. For g > the network contains flows which are 
cyclic both from the electric and the magnetic points of view, given by C(G) D C(G). Pop- 
ulating such a circular flow diminishes the ability of the network to respond to differences 
in potentials because it takes away from the pool of propagating charges. Hence the con- 
tribution of the modes e C(G) fl C(G) to the effective action is opposite from that of the 
zero modes of the fundamental field t/jq, which correspond to free excess charges. As r,, if) l 
are zero modes, this effect is not countered by a corresponding rescaling of the Laplacian 
eigenvalues, as is the case for all other modes. As a result, Tj, ipQ entirely capture the effect 
of conformal transformations on the effective action. This provides a natural interpretation 
of the familiar result 

/ 9au 6 , oc mdD with D = di + 6 1} (39) 

hi 09^ 

where D acts on the space of 1-forms Q l (M). In this language, the need for a doubling of 
the degrees of freedom, which is ubiquitous in geometric discretization [2], is understood as 



3 In particular, the theory presented herein is consistent even if the edge lengths defining the underlying 
discrete space do not satisfy the triangle inequality. 
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arising from the requirement that electric and magnetic degrees of freedom be treated in a 
duality-symmetric way. Geometric discretization without doubling treats the two kinds of 
degrees of freedom asymmetrically, viz. the paragraph below equation (l2~Ij) . which in itself is 
not inconsistent. Invariance under the electric-magnetic duality becomes necessary, however, 
after the introduction of the metric and the Hodge star. It is satisfying to have geometric 
discretization provide an intuitive understanding of the trace anomaly in both physical and 
algebro-geometric terms. 

The route toward discretization followed in Sec. [3] may also provide attractive settings 
for modelling emergent locality. The discrete counterpart of curvature proposed in f l37|) is 
not local in that it depends on the detailed form of the zero modes, which is determined by 
the full structure of the triangulation. This is again to be contrasted with the approach of 
Ko and Rocek [6], who exhibited the effective action fixed by demanding agreement with the 
conical (and thereby local) notion of curvature. A study of the discrepancies between the 
two results may reveal some lessons about the emergence of locality. 
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